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The main inference problems related to the Rayleigh distribution are the estimation of 
its parameter and the test of the hypothesis that a given set of observations is from such a 
distribution. It is shown that (in case of radio signals) the most efficient estimate of the 
parameter is obtained using the sample mean power. Complications may arise when data 
are missing or are autocorrelated. Methods are given to deal with such complications also. 



1. Introduction 

In a previous paper [Siddiqui, 1962] some problems 
of estimation and testing of hypotheses were dis- 
cussed in connection with the Rayleigh distribution. 
In the present survey article, for the sake of complete- 
ness, some parts of the previous paper will be re- 
peated. However, an attempt will be made to 
condense the material already covered. 

1.1. Notation 

The function P will be used generically for any 
distribution function and p for any density function. 
Thus P(x) and P(y), in general, will not be the same 
functions. If X is a random variable EX will 
stand for the expected value of X and var X for the 
variance "of X. Thus 



P(x)=Pr(X<x),EX=r 



xdP(x)> 



vnrX=EX 2 -(EX) 2 , 

and if P(x) is differentiable, P'(x)=p(x). 

2. General Properties of Rayleigh 
Distributions 

When an incoming plane wave (electromagnetic, 
sound, or some other kind) passes through a scat- 
tering medium, which at each instant of time may 
be considered as a sample from an ensemble of media 

with stable statistical characteristics, the outgoing 

-> 
wave R = (R cos A, R sin A) may be thought of as 

composed of many independent random wavelets 

-> 

R i =(R i cos A if Ri sin A t ), i=l, 2, . . . , n. The 

problem in its most generality will impose no re- 

strictions on the distributions of R t . However, under 
certain special circumstances, it may be reasonable 
to assume the following: 
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(1) Complete incoherence, which means that each 
of At is distributed uniformly over the interval 
(0, 2ir). 

(2) Independence of phase and amplitude, i.e., 
Ri and A t are independent random variables. 

(3) Absence of a dominant vector, defined by the 
condition of the Central Limit Theorem [Cramer, 
1951, pp. 215-216, 185-186]. If the first two condi- 
tions are satisfied, this third condition will be 
satisfied if: (E7%R a i) 1/d (EXRf)- m ->0 as n-*». 
Thus, in particular, if Ri have identical distributions 
or if Ri = r, a constant for all /', then (3) is satisfied. 

The condition (3) will not be satisfied if, say, a 
dominant vector is present. Under such circum- 
stances the Nakagami-Rice distribution applies 
[Nakagami, 1940; Rice, 1944, and 1945|. 

If the above assumptions hold R cos A and R sin 
A are asymptotically independent normal variates, 
each with mean zero and variance (1/2) o- 2 , where 
The distribution of (R, A) is approxi- 
iven by 



o*=ER 

matclv 



2r 



da 



dP(r,a)=^e- rZ/ °(lr^>0<r<oo j 0< a <:27r, 

and R and A are asymptotically independent. We, 
therefore, consider a random variable R, which has 
the probability density and distribution function, 
respectively 



p(r)=2<r- 



-^ 2 ,P(r)-l- 



2 ,0<r<co. 



This distribution is known as the Rayleigh am- 
plitude distribution. We immediately have p(Q) = 
P(od)=0 i and p(r) has a unique maximum at the 
mode r = 2~ l/2 a=0. 707a. The quantile function , r{p), 
that is, the inverse function of P(r) is given by 



r(p) = a[-\n(l-p)YI 2 , 0<p<l, 



(2.2) 



and gives the 100 p percentile, r(p), of the distribu- 
tion (2.1). Thus if [ — In (1— p)] 1/2 is plotted hori- 
zontally against r(p) vertically we obtain a straight 
line through the origin with slope a. The values 
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of r(p)/a are given in table 1 column 2, in particular , 
the median 



r(0.5)=cr(ln2) 1/2 =0.832(r 



For F>-2 
ER k 



- 2 --'f 



/ff2 dr=a k T(l+k/2). (2.3) 



The mean, ER, and the variance, var R, are 

#i?=(l/2)7r 1/2 cr=0.866(r, var R=(l-7r/4)a 2 =0.2U6o- 2 . 

Note that £7?>r(0.5)>r , which indicates that the 
right-hand tail of the probability density function is 
much longer than the left-hand one (positive-skew- 

ness) . 

Set Z=R 2 , y = a 2 . If R is identified as the ampli- 
tude of a wave, then Z is its power and y = EZ the 
mean power. The probability density and dis- 
tribution functions of Z are, respectively, 



p(z)=y- 1 e- z/y ,P(z) = l-e- z/y ,0<z<oo. 



(2.4) 



This distribution is known as Rayleigh power distri- 
bution or exponential distribution. The density func- 
tion has its maximum at z = and decreases mono- 
tonically to zero as z increases to infinity. If z(p) 
is the quantile function, then z(p) = r 2 (p), i.e., 



z(p) = -y\n (l-p),0<p<l. 



(2.5) 



If —In (1— p) is plotted horizontally against z(p) 
vertically we obtain a straight line through the 
origin with slope y. The values of z(p)/y are given 
in table 1 column 3. If k^> — 1 



E(Z/y) k = rz h e- z dz = T(k+l). (2. 
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Thus EZ=y, var Z=y 2 . Differentiating (2.6) j 
times with respect to k and then setting k = 0, we 
obtain 

#[ln(Z/7)]'=J (ln^-^=r^(l),j=0,l,... . 

Thus 

# (In Z) =ln 7+^(1) =lny--0.577 

var (In Z)=^'(l) = 7r76= 1.64493, 

where \f/(x) = T' (x)/T(x) . It is customary to measure 
power in decibels. Thus, writing (?=10 log 10 
Z=20 logic R, we have 

EQ=(10 log 10 e) E(ln Z) -4.343 (In 7-0.577) 



var #=(10 log 10 e) 2 var (In Z)=-31.0 



We note that the variance and also other central 
moments of Q are independent of the parameter y 
This is not unexpected since the distribution of 
(Z/y), hence that of In Z— 1117, is free of parameters. 

The quantile function, q(p), of Q is easily evaluated 
without explicitly writing down the distribution 
function of Q. In fact, the events "Q< q(p)" and 
"Z<z(P')" have the same probability if q(p) = lO 
logio z(jp). 



Thus 



q(p) = 10 log 10 7+10 log 10 (2(2?)/7). 



Hence q(p)—q(0.5) is free of the parameter 7. This 
quantity, which measures the distance of 1 the 100 
p percentile from the median, is tabulatedin table 
1 column 4. 

Table 1. Percentiles of R, Z, and Q 



(1) 


(2) 


(3) 


(4) 


lOOp 


r(p)/a 


z(p)/r 


q(p)-g(0.5) 


0.1 


0. 0316 


0. 00100 


-28.4 


1 


.100 


.0101 


-18.4 


5 


.227 


.0513 


-11.3 


10 


.324 


.105 


-8.18 


25 


.537 


.288 


-3.82 


50 


.832 


.693 


0.00 


75 


1.18 


1.39 


3.01 


90 


1.52 


2.30 


5.21 


95 


1.73 


3.00 


6.36 


99 


2.14 


4.61 


8.22 


99.9 


2.63 


6.91 


9.98 



A comparison of (2.4) with (2.1), or (2.5) with 
(2.2) indicates that the properties of the Rayleigh 
power distribution are much neater arid better 
suited for statistical estimation and for testing of 
statistical hypotheses. In fact, if we even start 
with (2.1), we find that the maximum likelihood 
(also the unbiased sufficient) estimator of a 2 from 
random sample R h . . ., R N is A^ 1 S?i??=iV- 1 
2 J iZ 2 , which is readily understood as (sample) 
"mean power" rather than "mean of amplitude 
square." Similarly, the reparametrization, with 7 
replacing a 2 , simplifies our estimates. From here 
onwards, therefore, we will mainly study the Rayleigh 
power distribution as given in (2.4). 

3. Estimation 

3.1. Estimation From a Random Sample 

The only essential parameter to be estimated in a 
Rayleigh power distribution is the mean 7. All other 
characteristics of the distribution such as the distri- 
bution function, momemts and percentiles are func- 
tions of 7. Let Z], . . ., Z N denote a complete 
sample of A^ independent observations from the dis- 
tribution. The likelihood function (the joint prob- 
ability density of Zi, . . ., Z N considered as a 
function of 7) is given by 

L(y)=y-» exp (-2Z t h). 



(2.7) I This shows that 2Z, a sufficient statistic for y 
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Since E2Z t =Ny, the unbiased 'sufficient (also the 
maximum likelihood) estimator of 7 is the sample 
mean, 

c=N- l VZ& (3.1) 

The variance of c is, of course, var c=y 2 /N, and if c' 
is any other unbiased estimator of 7 then 
var c<var c'. Since 2Z t /y is a x 2 variate with 2 deg 
of freedom and Z u . . ., Z N are independent, 

2Nc/y=22Z { /y 

is a x 2 variate with 2N deg of freedom. Thus 

X=Nc/y (3.2) 

has the distribution 

dP(x) = [T(N)]- 1 e- x ^- 1 dx } 0< x<°o, (3.3) 

from which the distribution of c is easily derived by 
substitution from (3.2). We note that if £> — TV 

EX t = [T(N)]- i rx t+N - 1 e- x (lx=T(t+N)/r(N). 

(3.4) 

Since c is the sufficient statistic for the family of 
Rayleigh distributions {P{z)\ 0<7<°° }, all other 
characteristics of the distribution should be estimated 
in terms of c. For example 

P(z) = l-e- z/c , 0<2<oo, (3.5) 

is the maximum likelihood estimator of P(z) and 

£(p)=-cln(l— p), 0<^<oo, (3.6) 

that of z(p). In general, if a function f(y) has an 
unbiased estimator, it also has an unbiased estimator 
based on the sufficient statistic c, and this estimator- 
is the most efficient estimator off(y). Also let h(x) 
be a function which in some neighborhood of the 
point x=y is continuous and has continuous deriva- 
tives of the first and second order. Then, if 
\h(x)\<Ce Kx for all x>0 where K is some constant 
independent of x, it can be shown that, for suffi- 
ciently large N, 



Eh(c)=h(y)+0(N- 1 ) 
var h(c)=N- l y 2 h'\y)+0(N-* /2 ). 



(3.7) 



For example if we wish to estimate the mean 
amplitude ER=(iry) l/2 /2, we consider m=(Tc) 1/2 /2 
and from (3.7) 



Em=(Try) 1/2 /2 + 0(N- 1 ),\xr m-- 



Try 



im 



-0(JV" 8/2 ). 



However, in this case, exact evaluations of Ec h 



and var c 1/2 are possible. In fact, from (3.2) and 
(3.3) 



2 _r(iV+i/2) 2 
^ c ~7v 1/2 r(JV) 7 ' 



w9 r T>(N + l/2) l 
~~L NT 2 (N) J 7 ' 



Thus an unbiased estimate of Eli and its variance 
are 

(*N)>*T(N) ._»/ Nl*(N) 

" 2T(N+l/2) 



var m" 



-/JVr ! 0V L _A 

!vT 2 (AT+l/2) J 7 ' 



(3.8) 

For small N these exact expressions can be used. 
For large N, 



m* 



= l/2(l+^)(7rc)'«+0(iV- 2 ), 



var m 



7T7 

; 16iV 



+0(iV- 2 ), 



which is an improvement over the biased estimate 
m._ On the other hand, let R=N- l ^R u then 
ER=ER, 

var n= w 7 * 

Thus the relative asymptotic efficiency, e(R), of R is 



7? N var mr 



= 0.91, 



' (jB) ~varl 4(4-tt) 

i.e., the sample mean R is a less efficient estimator 
of the population mean ER than the estimator m* 
based on the sufficient statistic c. 

If the measurements are in decibels, that is, if our 
sample consists of Q i =\Q log 10 Z t , i=l, 2, . . . ,N, 
then we note from (2.7) that an estimation of EQ, 
which involves In y, is required. We can directly 
evaluate 2?(ln c) and var (In c) by differentiating 
(3.4) once and twice with respect to J and then 
setting £=0. We note that 

E{\nXy=T«\N)lT{N), j=l, 2, . . ., 

^(iV)=r / (iV)/r(JV)=*(i)+£ 1 r- 1 , 

r = l 

+'(N)=ib (N+r)-°. 

Thus, since In c=ln X+ln y— In A 7 ', 

£" (In c)=ln 7-ln A^+^(vV)=ln 7 + 0(A^-') 

var (In c)=var (In X)=f'(N). 

Thus, to order JV _1 , the unbiased estimator of EQ 
based on c is 



#*=4.343 anc-0.577),varQ*=(4.343)V'(JV). (3-9) 
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On the other hand, if Q=N- 1 SQ t! EQ=EQ and 
7; (4.343tt) 2 

Thus the efficiency e N (Q) of Q relative to Q* is 
/7 ^_ var Q* _W . 



e N (Q)* 



var Q 



J2(N+r)- 



For N=l, e N (Q) = l, otherwise e N (Q)<Cl. Using 
the Euler-Maclaurin formula 

S M = (j(x)dx-i/2[f(k) -y(o) ] 

r-0 J° 



+j^L/ a) (*)-/" (0)]-... 



(3.10) 



to sum the series in 6^(0 we obtain, for large 
2V(>10), 

n 



Thus the efficiency of Q is approximately 61 per- 
cent. 

We have thus established the fact that when a 
complete random sample is available, one should 
invariably use estimates of population character- 
istics using appropriate functions of the sufficient 
statistic c. The use of the other estimators men- 
tioned above involves considerable loss of efficiency. 

3.2. Confidence Limits 

We noted that 2Nc/y is a x 2 variate with 2N deg 
of freedom. To set up confidence limits on y with 
a confidence coefficient 1 — a we determine from x 2 
tables two numbers a and b, corresponding to 2N 
deg of freedom, such that 

Pr( x 2 <a) = a/2, Pr{ x 2 <b) = l-al2. 
Then a<2Nc/y<b has probability I — a and 

2Nc/b<y<2Nc/a (3.11) 

are the desired confidence limits. If h is a monotone 
function from (0, oo), then h(2Nc/b) and h(2Nc/a) 
will be 100(1 — a) percent confidence limits for h(y). 
Thus, for example, EQ and P(z) are monotone func- 
tions of 7, and confidence limits for them can be 
easily written down. 

If iV>15, x=(4Nc/y) 1/2 -(4N-iy /2 is approxi- 
mately a standard normal variate [Cramer, 1951, 
p. 251]. If x a is the number such that 



& 



''-'°/.: 



e~ (1/2)x2 dx=l- 



then 4Nc[(4:N-iy /2 -±x a ]- 2 are 100(1 -a) percent 



confidence limits for y. If iV>100, these limits are 
approximately N 1/2 c(N 1/2 ±x a )~ l . 

3.3. Estimation From Order Statistics 

An estimate of y can be made with reasonably high 
efficiency from a few order statistics. Let Z l7 . . ., 
Z N denote, as before, a sample of N independent 
observations on Z and let Yi<Y 2 < . . . <Y N be 
their ordered values, e.g., Fi=min (Z u . . ., Z N ), 
Y N =msiX (Zj, . . ., Z N ). The probability density 
function of Y k is given by 



p(y) z 



N\ 



exp [— (N— 



(k-l)l(N-k)\y 

t+l)2//7][l-exp (-yMf- 1 , 

The moment generating function of Y k is 



0<y<oo. 



e uy p(y)dy= II (N—k+j)(N-k+j+yu)-\ 



From $ k (u) we easily evaluate, with a*= (N— i+ l)" 1 , 
EY k =yJ}a t , var F^ZX (3.12) 

i=l i=l 

Using the Euler-Maclaurin formula (3.10) to approxi- 
mate these summations, I have shown elsewhere 
[Siddiqui, 1963] that the optimum unbiased estimator 
from a single order statistic corresponds to 

£^0.79681 (iV+1) -0.39841 + 1.16312(AT+1)" 1 

(3.13) 

where the value of k thus obtained is rounded off to 
the nearest integer. For this k, 



Ci=IV2>^0.6275F, 

i 

is an unbiased estimator of y with efficiency 
var c 



(3.14) 



e N (ci)-- 



var Ci 



^.65+CKiV- 1 ) 



Comparing the efficiencies e N (c\) and e N (Q), we 
notice that even a single order statistic, properly 
chosen, is more efficient than the mean of all the 
decibel values for estimating the parameter of a 
Rayleigh distribution. It may also be noted that 
the efficiency of the estimator based on the median, 
i.e., the order statistic corresponding to fc^(l/2) 
(2V+1), is only 48 percent. 

The asymptotically (JV-^°°) optimum unbiased 
estimator of y from two order statistics is given by 

c 2 =0.523F m +0.179F n , m^0.639(iV+l), 

n^0.927(N+l), (3.15) 

with asymptotic efficiency e(c 2 )=0.82. Sarhan, 
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Greenberg, and Ogawa [1963] have given asymptoti- 
cally optimum unbiased estimators of 7 from 1 to 15 
order statistics. A linear combination of as few as 
five properly chosen order statistics provides an 
estimator of 7 whose efficiency exceeds 94 percent. 

A big advantage of the estimators based on order 
statistics is that only a partial knowledge of the 
sample is required. They are especially good when 
low values of the sample are poorly recorded, not 
recorded or missing. Let us say that the lowest 60 
percent of the sample is missing. We can still 
construct the estimators d and c 2 , although the 
sample mean (even median) can not be determined. 
The optimum estimators from more than one order 
statistic cannot be constructed if, say the highest 
8 percent or more values in the sample are missing. 
If only the highest values are missing so that we are 
left with Tij . . . , Y m , m<N, the most efficient 
unbiased estimator of 7 from these order statistics is 



-m 



[S Y t +(N-m)Y m l (3.16 



6) 



and has efficiency e(c m ) = m/N. In fact, 2mc m /y is a 
X 2 variate [Epstein and Sobel 1954] witli 2m deg of 
freedom. Confidence limits based on c m can be 
constructed in exactly the same way as with c. 

3.4. Estimation From Correlated Observations 

In this section we consider Z=Z(t) to be a station- 
ary Rayleigh process and assume that for all t and s 

EZ(t)=y,EZ(t)Z(t+s)=y 2 (l+a(s)), (3.17) 

so that a(s) is the correlation function of the process. 
Note that a(0) = l,a(—s) =a(s). Let Z(t) be observed 
over an interval of time 0<t< T. From the analogy 
of the random sample estimate, we use the sample 
mean 

c T = T~ 1 f Z(t)dt (3. IS) 

Jo 

to estimate 7. c T is unbiased and has variance 

vsirc T =2y 2 T- 1 f (l-s/T)a(s)ds. (3.19) 

We note that if we choose to take a discrete set of 
equally spaced observations Z(h), Z(2h), . . ., 
Z(Nh), then the mean, c N , will have the variance 



y 2 o*y 2 ^v-i / s \ 



0) 



The exact distribution of c T or c N is unknown. How- 
ever, let N' be denned by the equation 



var c=y 2 /N' , i.e., N'=y 2 /v&r c, 



(3.21) 



where var c is either (3.19) or (3.20). Then 2N'c 7 /y 
or 2N'c N /y is approximately a x 2 type variate with 
2N' deg of freedom. 



Generally, the correlation function a(s) will be 
unknown and will have to be estimated from the 
sample. Let 



Then 



C(s) = (T-s)- 1 f 1 " Z(t)Z(t+s)dt. 
a(s) = C(s)/c 2 T -l 



is a consistent estimator of a(s). Usually, it is 
desirable to fit to a(s) some mathematically specified 
function such as exp ( — ijl\s\), exp ( — /z.v 2 ), ^(ju+s 2 ) -1 . 
Before choosing an approximating function, however, 
a more careful investigation of a(s) near s = is 
necessary. For example, if we fake exp ( — /x\s\) to 
represent a(s) then, since this function is not differ- 
entiable at 6' = 0, we will be compelled to conclude 
that the process Z(t) is not differentiable, and this 
may not be a very desirable state of affairs. 

For the sake of simplicity we take the unit of time 
such that r=l. Lei 

yi =z(^)-z(^y i=i,2,...,N, 

where A 7 is given different values, say, 25, 50, 100, 
etc. We calculate 



1 
It can be shown that 

plim S N =y 2 [a' (-0) -a' (+()) \ 

iV-»oo 

where a { — 0) is the left-hand and a'( + 0) is the 
right-hand derivative of a{s) at s=0. If a(s) is 
differentiable at zero then this limit is zero, other- 
wise not equal to zero. Thus the behavior of S N as 
N increases will tell us whether a(s) is dilYerentiable 
at s—0 or not. Now supposing a(s) is twice differ- 
entiable at s = with a f (0)^(), we must have 



plim S N =0j plim NS N 



2 / / 
-7 a 



(0). 



We already know, in this case, that a'(0)=0. If 
exp ( — fxs 2 ) seems to represent a(s), then an estimate 
of fi is obtained as m=NS N /(2c 2 T )' 

4. Tests of Statistical Hypotheses 
4.1. Tests for the Distribution Function 

Let Z b . . . , Z N denote independent observations 
on a nonnegative variate Z. We wish to test the 
statistical hypothesis that Z has Rayleigh power 
distribution. We employ the x 2 goodness-of-fit 
test for this purpose. We first calculate c = N~ 1 ?%Zi\ 
then the numbers x t =— c In (l—pt), p i =i/m i i=l, 
2, . . . , m— 1, where m>5. From (3.6) we see that 
Xt are 100^9. percentiles of the Rayleigh power distri- 
bution with y = c. The expected number of observa- 
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tions in each of the intervals (0, x{), (x b x 2 ), . . . , 
(x m - u oo) is N/m. Let/ l7 . . . , j m be the number of 
observed Zj% falling in these intervals. The statistic 



_^ m Vi-Nlm) 2 _m 



= ST 



N/m 



~mfi-N (4.1) 



is asymptotically a x 2 variate with m — 2 deg of 
freedom under the hypothesis that Z is a Rayleigh 
variate. We preassign some critical probability 
level a( = 0.05 or 0.01, say) and from x 2 tables de- 
termine the number x« such that Pr(% 2 > x«) =<*• 
If the observed x 2< Cx«, we accept the hypothesis 
that Z is a Rayleigh variate, otherwise we reject it. 
If the sample is truncated at one or both ends, we 
use order statistics to estimate y and proceed as above 
changing the first and last intervals, if necessary, to 
read (0, x) and (y, 00) where £ is a number slightly 
less than the lowest recorded observation and y is a 
number slightly greater than the highest recorded 
observation. The expected number of observations 
in (0, x) is N{l—e- x,c ) and in (y, 00) is Ne~ y/C . We 
calculate 

v 2 X-v \J0 Je) 

X_Zj /. 

where f Q = observed, j e = expected frequency. 

In case of correlated observations it is necessary to 
multiply the value of x 2 thus obtained by N" /N 
where N" is the equivalent number of independent 

observations. Thus j&=(N"/N)2 (fo ~^ e) should 

Je 

be considered as an approximate x 2 variate with m— 2 
degrees of freedom. It is difficult to determine N ff . 
However, if a(s) is the correlation function of Z(t) 
process and the sample consists oiZQt), . . . , Z(Nh), 

N"IN^± a a»(jfc)]' =[l+2 Sr a 2 (jh)J\ 

appears to be a reasonably good approximation to 
N"IN. 

4.2. Comparison of Two Samples 

Sometimes we wish to test the hypothesis that 
two independent samples, which are known to be 
from Rayleigh power distributions, are from the 



same distributions. Let C\ and c 2 be the means of the 
samples of sizes Ni and N 2 from Rayleigh distribu- 
tions with means y h and y 2 , respectively. Then, 
under the hypothesis y x =y 2} 

F(2N u 2N 2 )=c 1 /c, 

is a Fisher-Snedecor F variate with the indicated 
degrees of freedom. We assume that Ci>c 2 \ if not 
we simply invert the ratio and interchange the 
degrees of freedom, preassign a significance level a 
and test for the significance of the calculated F by 
comparing it with the upper 100a:/2 percentage 
point of the F distribution. If the observations in 
each sample are autocorrelated, we modify the 
degrees of freedom from N to N / according to (3.21). 
If the i^-test indicates acceptance of the hypothesis 
71=72=7, the common mean, y, is estimated by 
c=N~ 1 (N 1 c 1 +N,c 2 ),N=N 1 +N 2 . 
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